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At  the  request  of  the  U.S.  Amy  Missile  Research  and  Development 
Comnand  and  the  Harry  Diamond  Laboratories  the  Ballistic  Research 
Laboratory  has  undertaken  the  development  of  an  appropriate  methodology 
for  assessing  the  vulnerability  of  solid  rocket  motors  to  lateral  air 
blast  loading,  especially  that  produced  by  nuclear  weapons  detonation. 

It  is  desired  that  this  methodology  be  compatible  with  analysis 
procedures  currently  used  for  design  or  evaluation  of  rocket  motors, 
specifically: 

(1)  a finite  element  structural  analysis  code^,  AMG074,  for 
operational  loads  analysis, 

2 

(2)  a finite  element  stress  wave  analysis  code  , AMG063,  for 
treating  x-ray-induced  stress  waves, 

(3)  a photon  energy  deposition  code  , KNISH,  for  predicting 
energy  absorption  and  subsequent  stress  wave  loadings,  and 

4 

(4)  a computer  program  for  combined  operational  and  nuclear 
effects  loads  which  employs  the  three  preceding  codes. 

To  achieve  this  objective,  a contract  has  been  assigned  to  personnel 
familiar  with  the  cited  analysis  procedures  for  the  development  of  a 
two-dimensional  plane  strain  analysis  of  the  motor  cross  section, 
including  the  viscoelastic  propellant  grain,  elastoplastic  deformation 
of  the  motor  case,  effects  of  internal  pressurization,  operational  loads 
plus  side-on  blast  loading  subsuming  diffraction  effects. 


J.  J.  Briabanej  "Advances  in  Stress  Analysis  of  Solid  Propellant 
Rocket  Grains Rohm  and  Haas  Company  Technical  Report  S-268, 

Contract  DAAH01-7Q-C-01S2,  September  1970. 

2 

R.  S.  Dunhamj  "Dynamic  Analysis  of  One  Dimensional  Thermorheologically 
Simple  Viscoelastic  Solids  With  Nonlinear  Heat  Conduction  Analysis" y 
U.S.  Army  Missile  Command  Report  No.  RK-TR-70-1Z,  July  1970. 

7 

J.  Case  et  aly  "The  KNISH/AMGO  Photon  Energy  Deposition  and 
Viscoelastic  Material  Response  Computer  Program"y  Physics 
International  Companyy  PITR-307-Zy  August  1972. 

4 

T.  L.  Cost  and  G.  E.  WeekSy  "Automated  Evaluation  of  Interceptor 
Rocket  Motor  Designs  Under  Combined  Operational  and  Nuclear  Effects 
Load8"y  University  of  Alabama  BER  Report  No.  176-97 y July  1974. 
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However,  it  was  felt  desirable  to  explore  the  utility  of  the  BRL 

version  of  the  PETROS  3.5  code^  for  predicting  the  response  of  rocket 
motors.  This  computer  program,  which  is  based  on  a nonlinear  finite 
difference  analysis,  has  been  extensively  employed  to  treat  blast- 
induced  response  of  thin  elastoplastic  shell  structures.  By  contrast 
with  the  finite  element  analysis  model  to  be  produced  by  a contractor, 
the  PETROS  3.5  model  is  valid  for  large  deflections  and  can  be  employed 
to  treat  an  entire  rocket  motor  case  rather  than  just  a sectional  cut. 
However,  it  has  no  provision  for  inclusion  of  the  propellant  grain  in 
the  calculations  and  rather  extensive  modifications  to  the  code  would 
be  fequired  to  incorporate  even  plane  strain  deformation  of  the 
propellant.  Consequently,  in  this  report  the  PETROS  3.5  code  will  only 
be  used  to  treat  the  plane  strain  response  of  a bare  cylinder  (motor 
case),  which  may  in  the  future  be  compared  with  the  contractor's  finite 
element  analysis,  and  the  cylinder  plus  an  approximate  modeling  of  the 
mass  effect  of  the  propellant. 

The  representation  employed  for  the  blast  loading  includes  both  the 
diffraction  and  the  drag  phases  of  loading  since  both  phases  contribute 
to  the  response  to  nuclear  blast.  Nevertheless,  the  methodology  is 
equally  applicable  for  predicting  response  to  blast  from  conventional 
weapons . 

In  addition  to  results  of  numerical  analysis  obtained  with  the 
PETROS  3.5  code,  this  report  contains  analytical  solutions  for  the 
effect  of  pressurization  upon  the  vibration  characteristics  of  elastic 
cylinders  and  for  the  response  of  such  cylinders  to  transverse  blast 
loading.  There  is  also  a discussion  of  conditions  under  which  solid 
propellant  motors  are  most  susceptible  to  blast-induced  damage. 


II.  PLANE  STRAIN  CYLINDER  RESPONSE  CALCULATIONS 

When  resorting  to  numerical  analysis  one  must  naturally  assign 
values  to  all  parameters  appearing  in  the  formulation.  The  values 
selected  below  are  not  identified  with  any  existing  system  but  do  fall 
within  ranges  of  practical  interest. 

A.  Rocket  Motor  Parameters 


A representative  solid  propellant  rocket  motor  cross  section  (see 
Figure  1)  was  selected  for  exercise  of  the  PETROS  3.5  computer  program 
in  its  plane  strain  option.  In  addition  to  the  geometric  data  shown 


^S.  D.  Pirotin^  B.  A.  Berg,  and  E.  A.  WitmeVy  "PETROS  3.5:  New 
Developments  and  Program  Manual  for  the  Finite- Difference  Calculation 
of  Large  Elastic-Plastic  Transient  Deformations  of  Multilayer  Variable- 
Thickness  Shells",  U.S.  Army  Ballistic  Research  Laboratories  Contract 
Report  No.  212,  February  1975.  (AD  #A007215) 


Figure  1.  Rocket  Motor  Cross  Section 

on  this  Figure,  the  material  property  parameters  listed  in  Table  1 were 
selected.  The  steel  motor  case  was  assumed  to  have  the  linear  strain 
hardening  uniaxial  stress-strain  characteristics  shown  in  Figure  2 and 
to  exhibit  no  strain-rate  effects.  Bi-axial  plastic  flow  was  analyzed 
using  the  von  Mises  yield  criterion  and  the  associated  flow  rule  in 

accordance  with  the  mechanical  sublayer  model  . Although  the  propellant 
is  being  modeled  as  a viscoelastic  material  in  the  cited  finite  element 


L.  Marino t J.  W.  Leeahj  and  E.  A.  Witmerj  "PETROS  2:  A New  Finite- 
Difference  Method  and  Program  for  the  Calculation  of  Large  Elastic- 
Plaetic  Dynamically -Induced  Deformations  of  General  Thin  Shells", 
U.S.  Army  Ballistic  Research  Laboratories  Contract  Report  No.  12, 
December  1969.  (AD  #708773) 
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Table  I . Material  Properties 


Motor  Case  (Steel) 


Young's  modulus 

E « 200  GPa  * 
c 

29  X 10®  psi 

Poisson's  ratio 

V = 0.3 
c 

Yield  stress 

0 = 483  MPa  = 

0 

70,000  psi 

Mass  density 

= 7850  kg/m^ 

= 0.000735  lb  sec^/in'* 

Propellant  (Viscoelastic  Solid) 

Young's  modulus 

E = 689  MPa  = 

p 

100,000  psi 

Poisson's  ratio 

V = 0.472 

p 

Mass  density 

Pp  = 1660  kg/m^ 

= 0.000155  lb  sec2/in‘' 

Figure  2.  Postulated  Uniaxial  Stress-Strain  Curve  for 
Motor  Case  Material 

analysis,  the  lumped  mass  treatment  of  the  propellant  in  the  PETROS  3.5 
calculations  required  use  of  only  the  propellant  density. 
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B.  Blast  Loading  Model 


The  blast  loading  is  introduced  as  a plane  shock  wave  having  an 
exponentially  decaying  tail,  moving  down  from  above  as  shown  in  Figure  1. 


Although  surface  overpressures  p (e,t)  obtained  experiment ally, or  by 


refined  hydrodynamic  calculations  can  be  readily  introduced  into  the 
PETROS  3.5  code,  it  is  convenient  to  employ  the  following  functional 
representation  for  the  blast  loading  which  includes  the  essential  aspects 
of  wave  reflection,  diffraction,  and  post-envelopment  drag  loading: 


P (0.t)  = I [p  COS0  + p (l-cose)]e’^^^'^a^  for  -90®  < 9 < 90® 
c I r s 

for  90®  < |e|  < 180' 


P e 


where 


p = ambient  pressure 


p = side-on  overpressure  of  incident  shock 


p^  = reflected  peak  overpressure 


a = sonic  velocity  at  ambient  conditions 
o 


U = shock  front  velocity 


t = arrival  time  of  shock  front 
a 


R = outside  radius  of  cylinder 
X = decay  coefficient 


and 


Pr  = 2p, 


VPo  ^ 


4P. 


^ R(1  - cosB) 

" U 


for  t < t 


for  t > t 


(1) 


(2) 


(31 


(4) 
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While  response  predictions  have  been  made  for  several  values  of 
incident  overpressure,  the  results  which  follow  have  been  obtained  for 
one  level  of  blast  loading,  the  parameters  for  which  are  listed  in 
Table  II. 


Table  II.  Blast  Wave  Parameters 


Parameter 

SI  Units 

English  Units 

a 

0 

340  m/s 

13400  in/sec 

Po 

101  kPa 

14.7  psi 

Ps 

79.3  kPa 

11.5  psi 

U 

440  m/s 

17300  in/sec 

X 

4.35  s-1 

4.35  sec'^ 

Pr 

207  kPa 

30.0  psi 

The  value  of  the  decay  coefficient  X was  estimated  from  blast  data 

7 

curves  for  a one  kiloton  nuclear  weapon. 

C.  Internal  Loading 

In  addition  to  response  predictions  using  the  foregoing  external 
loading,  which  would  correspond  to  the  pre-ignition  response  for  the 
rocket  motor  stage  under  consideration,  calculations  have  also  been 
performed  for  cases  where  the  motor  had  an  additional  internal 
pressurization  p^  (produced  by  propellant  combustion)  prior  to  this 

blast  loading.  For  the  most  part  the  value  p^  = 6.89  MPa  (1000  psi) 

has  been  used  for  this  internal  pressure,  although  a limited  study  of 
the  effect  of  varying  this  parameter  will  also  be  reported. 


III.  FINITE  DIFFERENCE  MODELING 
A.  Computational  Grid 

The  finite  difference  grid  employed  for  all  the  PETROS  3.5 
calculations  is  shown  in  Figure  3.  Since  both  the  responding  structure 
and  the  applied  loads  are  symmetric  with  respect  to  the  vertical  plane, 


7 

S.  Glaaatone  (Editor) , "The  Effects  of  Nuclear  Weapons",  United  States 
Atomic  Energy  Commission,  April  1962, 


CYCLE  0 

Omt 


Figure  3.  Finite  Difference  Grid 

it  is  only  necessary  to  model  the  response  of  one-half  of  the  cylinder. 
Further,  because  there  is  no  longitudinal  variation  of  any  quantities 
involved  in  a plane  strain  analysis  there  are  only  two  independent 
variables,  the  angles  9 and  time,  resulting  in  quite  economical  co nputer 
runs. 

B.  Propellant  Mass  Effect 

In  the  finite  difference  equations  of  motion  for  the  mesh  point  at 

6 (m  = integer)  the  mass  is  that  of  the  shell  (motor  case)  lying 
m 

between  the  radial  planes  at  angles  (6  , + 6 )/2  and  {0_  + 6„^i)/2. 

m-l  m m m+i 

When  it  was  desired  to  take  account  of  the  mass  of  the  propellant  in 
the  response  prediction  the  mass  of  propellant  lying  between  the  same 
two  radial  planes  was  added  to  the  motor  case  mass  at  0^.  Therefore, 

the  effective  mass  will  vary  as  a function  of  the  discrete  angles  0^ 

owing  to  the  slots  in  the  star  pattern  of  the  propellant  grain  (see 
Figure  1) . It  is  recognized  that  this  procedure  provides  a rather 
crude  approximation  to  the  propellant  mass  effect  since  many  two- 
dimensional  response  modes  are  thereby  inhibited.  Nevertheless  it 
appears  to  be  the  most  rational  procedure  which  does  not  entail 
extensive  modification  of  the  shell  response  code.  The  resistance  of 
the  propellant  grain  to  deformation  is  also  neglected  in  this  modeling, 
an  approximation  which  may  not  be  too  serious  owing  to  the  relatively 


low  strength  of  the  propellant  in  comparison  to  that  of  the  steel  motor 
case  (see  Table  I)  and  the  weakening  effect  of  the  slots  in  the 
propellant  grain. 

C.  Rigid  Body  Motion 

There  is  another  modeling  consideration  which  arises  as  a consequence 
of  the  assumption  that  the  motor  case  is  in  a state  of  plane  strain.  In 
this  stress  state  there  can  be  no  resultant  force  acting  tangential  to 
a cross  section  of  the  motor.  Consequently,  when  a non-self-equilibrating 
external  loading  (such  as  results  from  side-on  blast)  is  applied  to  the 
cylinder  a lateral  acceleration  results.  In  order  to  appreciate  the 
extent  of  deformation  the  cylinder  experiences  it  is  necessary  to 
subtract  out  this  "rigid  body"  component  of  the  motion.  This  was 
accomplished  in  the  PETROS  3.5  calculations  by  recomputing  the  location 
of  the  mass  center  at  each  time  step  and  subtracting  this  quantity  from 
the  position  vector  of  each  mesh  point  before  displaying  graphical 
results. 


IV.  NUMERICAL  ANALYSIS  RESULTS 
A.  Unpressurized  Cylinders 

In  Figure  4 the  undeformed  cylinder  cross  section  of  the  bare  motor 
case  is  compared  with  that  at  30  ms,  which  is  approximately  the  time  of 
maximum  deflection.  It  should  be  emphasized  that  in  this  figure  the 


INITIAL 

CONFIGURATION 


CYCLE  1250 
30  ms 

Figure  4.  Deflection  Pattern  for  Unpressurized  Bare  Cylinder 
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deflections  are  plotted  to  the  sane  scale  as  the  initial  configuration; 
i.e.,  the  deflections  are  large,  entailing  both  nonlinear  geometric 
effects  and  elastoplastic  material  behavior.  The  deflected  cross 
section  corresponds  principally  to  the  lowest  frequency  flexural  mode 
for  a ring  or  cylinder  (the  n * 2,  s * 1 mode,  see  Chapter  V.A,).  The 
deflection  responses  of  a point  on  the  crown  line  (6  = 0“)  of  cylinders 
having  no  internal  pressurization  are  shown  in  Figure  5,  for  both  the 
bare  motor  case  and  the  motor  configuration  of  Figure  1 (where  the 
effect  of  propellant  mass  has  been  included  as  described  above).  As 
expected,  the  effect  of  the  propellant  is  to  reduce  both  the  amplitude 
and  frequency  of  the  response. 


Figure  5.  Deflections  at  6 = 0®  for  No  Internal  Pressurization 

The  predicted  circumferential  strains  at  the  crown  line  are 
displayed  in  Figure  6.  Since  the  curves  for  the  inner  and  outer 
surfaces  (of  the  motor  case)  are  nearly  symmetric  with  respect  to  the 
zero  strain  line  it  may  be  inferred  that  the  response  is  mostly  flexural, 
which  is  consistent  with  Figure  4.  Similar  strain  variations  were 
obtained  at  6 * 90®  and  180®.  At  192  ms,  by  which  time  all  plastic 
deformation  had  ceased,  the  external  load  was  removed  and  an  artificial 
damping  was  introduced  to  dissipate  rapidly  the  remaining  kinetic 
energy.  In  this  manner  it  was  possible  to  identify  the  values  of 
residual  deflection  and  residual  strains  indicated  on  Figures  5 and  6, 
respectively. 
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DEFLECTION  (inches) 


TIME  (ms) 

Figure  6.  Circumferential  Strains  at  6 = 0®  for 
No  Internal  Pressurization 

For  the  motor  case  containing  propellant  it  is  also  of  interest  to 
estimate  the  normal  stresses  acting  at  the  interface  between  motor  case 
and  propellant.  This  can  be  done  in  a manner  consistent  with  the 
lumped  mass  model  by  multiplying  the  element  of  propellant  mass 
assigned  to  a nodal  point  by  the  absolute  acceleration  of  the  nodal 
point  and  dividing  by  the  appropriate  interface  surface  area.  The 
acceleration  histories*  of  the  nodal  points  at  6 = 0®  and  0 = 180®  are 
shown  in  Figures  7 and  8,  respectively.  Tensile  interface  stresses  are 

The  presence  of  high  frequency  oscillations  in  these  plots  may  appear 
suspicious  in  view  of  the  raider  low  frequency  deflection  response  in 
Figure  S.  However^  relatively  low  amplitude  extensional  modes 
(the  s = 2 family  discussed  in  Chapter  V)  when  amplified  by  the  square 
of  their  large  circular  frequencies  (as  is  effectively  done  when 
calculating  accelerations)  can  become  dominant  in  the  manner  shown  in 

g 

Figures  7 and  8.  Love  remarks  that  such  modes  would  probably  be 
difficult  to  excite.  It  appears  that  diffracting  a shock  wave  around  a 
cylinder  is  quite  effective  in  exciting  these  high  frequency  modes. 

Q 

A.  E.  H.  Lovsj  "A  Treatise  on  the  Mathematical  Theory  of  Elasticity''^ 
Fourth  Edition^  Dover  Publications ^ 1944 ^ p.  454. 
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associated  with  positive  accelerations  at  6 ■ 0°  and  negative  accelera- 
tions at  6 ■ 180*.  The  maxinum  of  these  accelerations  is  -496  m/s^ 
(-19500  in/sec^)  corresponding  to  a tensile  stress  of  80.2  kPa  (11.6  psi) 
at  6 ■ 180*.  While  interface  stresses  of  this  magnitude  should  not  cause 
debonding  of  the  propellant,  conclusions  regarding  bond  failure  should 
be  deferred  until  the  results  of  the  two-dimensional  finite  element 
analysis  become  available. 

B.  Effect  of  Internal  Pressure 


When  treating  cases  of  cylinders  with  internal  pressurization  (due 
to  propellant  combustion)  a quasi-static  PETROS  3.5  run  was  made  at 
each  pressurization  level  to  establish  the  initial  conditions  of 
deflection  and  stress  extant  prior  to  arrival  of  the  blast  wave.  These 
quasi-static  solutions,  which  are  plotted  in  negative  time  in  Figure  9, 

are  fully  damped  before  the  blast  loading  is  initiated  at  time  = 0.  The 

transient  responses  to  the  standard  blast  pulse  are  plotted  in  positive 
time  for  three  cases:  the  motor  bottle  with  and  without  propellant 
having  an  internal  pressure  of  6.89  MPa  (1000  psi),  which  is  understood 
to  be  a representative  operating  pressure  for  a rocket  motor,  and  the 
motor  without  propellant  with  half  the  above  internal  pressure  which  was 

p;  «6.89MPo  (1000  psi);  WITHOUT  PROPELLANT 

p.  . 6.89 MPo  (1000  psi  );  WITH  PROPELLANT 

p.  . 3.45 MPa  ( 500  psi ) ; WITHOUT  PROPELLANT 


DEFLECTION  (inches) 


I 


1 


run  to  study  the  effect  of  varying  the  internal  pressure.  One  sees  that 
increasing  the  internal  pressure  raises  the  frequency  of  oscillation 
and  decreases  the  amplitude  of  response  (provided  the  shell  remains 
elastic).  At  the  operating  pressure  adding  propellant  lowers  the 
frequency  of  response  and  has  little  effect  on  the  amplitude. 

Figure  10  presents  cross  section  views  of  cylinder  response  at  two 
times  for  the  case  represented  by  the  solid  curve  in  Figure  9.  In 
these  views  the  departure  from  the  initial  unpressurized  configuration 
has  been  magnified  by  a factor  of  one  hundred.  For  Cycle  0 we  see  the 
axisymmetric  extensional  deformation  pattern  which  exists  immediately 
before  the  arrival  of  the  blast  loading.  With  the  6.89  MPa  (1000  psi) 
internal  pressure  the  circumferential  membrane  stress  is  91  percent  of 
the  uniaxial  yield  stress  while  the  maximum  bending  stresses  are  less 
than  one  percent  of  this  quantity.  The  cross  section  for  Cycle  377, 
which  may  be  associated  with  the  second  inward  peak  displacement  shown 
in  Figure  9,  exhibits  an  essentially  inextensional  flexural  deformation 
superposed  upon  the  extensional  pattern  for  Cycle  0.  The  corresponding 
quasi-static  and  transient  strains  at  the  crown  line  are  depicted  in 
Figure  11.  Clearly,  the  blast-induced  strain  fluctuations  are  small  in 
comparison  to  the  pressurization-induced  strain.  In  fact,  in  spite  of 
the  large  pre-stressing,  none  of  the  problems  treated  in  Figure  9 were 
predicted  to  experience  permanent  deformation  of  the  motor  case  as  a 
result  of  the  blast  loading. 


INITIAL  UNPRESSURIZED  CYLINDER 


CYCLE  0 

0 mt 


CYCLE  377 
9.048  ms 


Figure  10.  Cross  Sections  of  Pressurized  Cylinder 
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Figure  11.  Circumferential  Strain  at  6 = 0®  for 

= 6.89  MPa  (1000  psi)  Without  Propellant 


The  stiffening  effect  of  internal  pressurizatio.T  is  demonstrated 
in  Figure  12,  which  shows  the  variation  of  maximum  deflection  (from  the 
pressurized  state)  of  the  bare  cylinder  as  a function  of  internal 
pressure  for  a constant  external  blast  loading.  The  black  dots  on  this 
figure  are  values  obtained  by  PETROS  3.5  calculations  and  the  curve  is 
a line  faired  through  these  values.  The  curve  is  dashed  beyond 
Pj^  = 8.48  MPa  (1230  psi)  to  indicate  that  a different  trend  may  be 

expected  after  the  internal  pressure  alone  is  sufficient  to  produce 
yielding.  The  great  magnification  of  response  in  the  neighborhood  of 
Pj^  = 0 is  of  considerable  practical  significance.  This  subject  will  be 

elaborated  upon  in  the  following  chapter. 


V.  ANALYTICAL  MODELING  AND  INTERPRETATIONS 

The  preceding  results,  which  were  obtained  by  numerical  solution 
of  the  nonlinear  partial  differential  equations  describing  finite 
amplitude  elastoplastic  response  of  thin  shells,  can  be  usefully 
interpreted  by  comparison  with  the  deductions  from  an  exact  analysis  of 
a linearized  formulation  for  the  same  problem.  Although  several 
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Figure  12.  Effect  of  Internal  Pressure  on  Peak  Response 
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MAXIMUM  DEFLECTION  ( inches) 


alternate  versions  of  the  latter  formulation  were  investigated  it 
appears  that  the  simplest  theoretical  development  which  contains  all 
the  desired  features  can  best  be  presented  as  a generalization  of 

M P •• 

Flugge's  equations  for  the  buckling  of  cylinders.  When  Flugge's 
formulation  is  particularized  to  the  plane  strain  case  and  augmented 
by  the  addition  of  inertia  and  load  terms  there  results: 


The  symbols  employed  in  these  equations  are  defined  in  the  List  of 
Symbols  and  the  sign  conventions  for  positive  values  of  displacements 
and  loads  are  indicated  in  Figure  13.  It  should  be  reiterated  that  this 
analysis  is  valid  only  for  elastic,  small  deflection  response. 


Figure  13.  Sign  Conventions  for  Displacements  and  Loads 

^W.  Flugget  ''Stresses  in  Shells"t  Second  Edition,  Springer-Verlag,  1973, 
p.  448. 
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If  desired,  one  can  eliminate  the  tangential  displacement  v and 
obtain  a single  differential  equation  for  the  radial  displacement  w. 
To  accomplish  this,  Equation  (5)  is  differentiated  with  respect  to  0 

3v 

and  Equation  (6)  is  solved  for 


(1  + q) 


i 

0 ieat^ 


(7) 


^ = -1— [-  w 

39  1 + q L 


+ 2 


w) 


e 3t^_ 


(8) 


Successive  differentiation  of  Equation  (8)  yields  the  terms  needed  for 
substitution  in  Equation  (7) when  this  is  accomplished  and  terms  are 
collected  one  obtains: 


a^w 

3^ 


(-!)«?•(■  ■?)«;• -S' 


6(1 

A.  Free  Vibrations 


Before  attempting  a solution  of  Equation  (9)  for  the  forced 
response  of  the  cylinder  it  is  expedient  to  consider  the  free  vibration 
characteristics  of  the  system  defined  by  the  preceding  equations. 
Setting  p = 0 in  Equation  (9)  and  substituting 


w 


cos  n0 


sin  oj  t 
n 


(10) 


demonstrates  that  the  function  cos  n0  is  an  eigenfunction*  and  leads  to 
the  frequency  equation: 


o^=||.|l.f(l-».2q)n2+k(l-n2)^±  [(l+n2)^+4(2+q)qn2+2k(l-n2)  -^k^Cl-n^)  ] 


(11) 


n takes  on  integral  values.  Sin  nd  is  also  an  eigenfunction  but,  since 
we  are  concerned  only  with  modes  symmetric  with  respect  to  9=0,  no  use 
will  be  made  of  this  function. 
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Equation  (H)  reveals  that  there  are  two  natural  frequencies 
associated  with  each  value  of  n.  Consequently  in  the  sequel  the 
frequencies  will  carry  two  subscripts;  e.g.,  os^  where  s = 1 or  2. 

Further,  it  should  not  be  assximed  that  the  v and  w components  of 
displacement  are  independent.  If  the  radial  component  of  the  nth  mode 
is  assumed  to  be 

w a cos  n0  (w  , sin  u ,t  + w - sin  w -t)  (12) 

n ^ n,l  n,l  n,2  n,2  '' 

where  w^  ^ and  w^  ^ ®re  arbitrary  amplitude  coefficients,  then  the 
tangential  displacements  in  the  same  mode  must  have  the  form 

/ 2 3 4\V1 

- <(n^+l)  ♦ 4(2+q)qn^  - 2k(n^-l)  + k^Cn^-l)  j Jw^  ^ sin 

+ j^(n2-l){l-k(n2-l)} 

+ |(n2+l)  + 4(2+q)qn^  - 2k(n2..1)  + k^Cn^-l)  | 2 *"n,2^| 

in  order  that  v , w satisfy  Equations  (7)  and  (8)  simultaneously. 

Since  either  of  the  amplitude  coefficients  w , or  w _ can  be  taken  to 

n,l  n,2 

be  zero,  it  is  clear  that  there  are  two  distinct  families  of  modes 
having  frequencies  ^ and  2*  lower  frequency  family  being 

identified  with  s = 1. 

Certain  useful  characterizations  of  these  two  families  of 
vibration  modes  can  be  achieved  by  substitution  of  their  individual 
displacement  components  into  the  expressions  for  circumferential 
membrane  strain 


= -fe  - A 

ayae  J 


and  for  circumferential  bending  moment 


M _ Eh! 

^0  12(l-v 


-v2)a2  ) 
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For  the  s « 1 modes  with  n > 2 the  membrane  strains  are  insignificantly 
small  while  the  bending  moments  are  appreciable.  Thus,  this  lower 
frequency  set  of  modes  can  be  identified  as  the  family  of  inextensional 
bending  modes  of  a cylinder.  The  special  cases  for  n >:  0 and  n s 1 
(for  s ■ 1)  are  not  vibratory  modes;  n ■ 0 corresponds  to  a static 
expansion  or  contraction  of  the  cylinder  while  n = 1 represents  a 
rigid  body  translation  of  the  cylinder  cross  section  along  the  6=0 
axis. 

The  s = 1 family  of  modes  is  extremely  sensitive  to  the  internal 
pressure  p^  as  illustrated  in  Figure  14  for  the  lower  natural  frequencies 

f , (=  <i)  ,/27r).  The  frequency  of  the  fundamental  inextensional  mode 

y X n ^ x 

f-  ^ is  seen  to  have  an  intercept  on  the  p.  axis  at  -0.214  MPa 

(-31.1  psi)  corresponding  to  the  value  of  uniformly  distributed 
external  pressure  for  which  static  buckling*  of  the  cylinder  would  occur 
in  the  cos  20  mode.  The  great  magnification  of  response  in  the 
neighborhood  of  p^  = 0 exhibited  in  Figure  12  can  be  explained  on  this 


Pj  (P»i) 

'200  0 200  400  600  800  1000  1200 


Pi  (MPo) 

Figure  14.  Effect  of  Internal  Pressure  on  Certain 
Inextensional  Mode  Frequencies 


See  Refevenoe  9,  pp,  459-60. 
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basis.  The  actual  blast  loading  can  be  decomposed  into  a transient 
uniformly  distributed  external  pressure  and  a non-uniform  transient 
loading.  The  latter  serves  to  deform  the  cylinder  from  the  initial 
circular  cross  section  (analogous  to  a large  initial  imperfection) 
while  the  uniform  coiq>onent  provides  a nonlinear  amplification  of 
response  if  its  magnitude  is  a significant  fraction  of  the  pressure 
required  for  static  buckling.  For  the  blast  loading  employed  in  the 
PETROS  3.5  calculations  the  temporal  variation  of  the  uniform  component 
is  shown  in  Figure  15.  This  pressure  component  reaches  a maximum  of 
0.119  MPa  (17.3  psi) , which  is  56%  of  the  static  buckling  pressure,  and 
decays  slowly  with  respect  to  the  period  of  the  fundamental  bending 
mode  (39.3  ms)  so  there  is  ample  time  for  a large  deflection  to  occur. 
The  response  for  p^^  = 0 was  also  enhanced  by  the  increased  flexibility 

associated  with  plastic  flow.  For  large  internal  pressurization  the 
circumferential  membrane  stresses  remain  tensile  in  spite  of  the  blast 
loading  so  there  is  no  near-buckling  effect.  In  fact,  consonant  with 
the  raising  of  the  frequency  spectrum  of  the  inextensional  modes  with 
increasing  p^^,  as  shown  in  Figure  14,  there  is  the  decrease  in  response 

amplitude  illustrated  in  Figure  12  (at  least  until  p^  alone  causes 

plastic  yielding). 
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Figure  15.  Variation  of  Mean  External  Pressure 
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UNIFORM  PRESSURE  COMPONENT  (psi) 


An  interesting  correlation  between  analytical  and  numerical 
predictions  is  shown  in  Figure  16.  The  circled  points  represent 
periods  j (“  2tt/uj2  j)  of  post -peak  elastic  oscillations  measured 

from  plots  of  PETROS  3.5  calculations  for  various  value •»  of  p^.  The 

solid  curve  is  the  analytical  prediction  of  variation  of  period  with 
internal  pressure  obtained  with  the  aid  of  Equation  (11).  The  agreement 
is  amazingly  good  for  the  two  larger  values  of  p.  and  the  deviations 

at  the  two  lower  internal  pressurizations  can  be  readily  explained.  At 
these  lower  values  of  p^  the  mean  value  of  the  external  blast  loading 

depicted  in  Figure  16  is  large  enough  relative  to  p^  to  produce  a small 

negative  shift  in  the  effective  value  of  the  abscissas  for  the  circled 
points  (which  would  produce  better  agreement  with  the  solid  curve). 

Pj  ( P*i) 

-200  0 200  400  600  800  1000  1200 


To  verify  this  deduction  the  following  tests  were  performed.  At 
= 0.345  MPa  (50  psi)  the  PETROS  3.5  code  was  run  to  obtain  static 

equilibrium,  then  instead  of  applying  the  blast  loading,  the  numerical 
solution  was  restarted  with  an  initial  velocity  distribution  corresponding 
to  the  (2,1)  mode  with  an  amplitude  such  as  to  produce  a response 
approximately  equal  to  that  produced  by  the  blast  loading.  The  post- 
peak vibratory  period  was  found  to  be  that  surrounded  by  the  small 
triangle.  At  p^  « 0 a PETROS  3.5  run  was  made  with  the  blast  loading 

scaled  down  to  a peak  value  of  6.89  kPb  (1  psi)  so  that  the  mean 
component  of  the  blast  loading  was  negligible  and  the  response  amplitude 
was  small.  This  led  to  the  free  vibration  period  enclosed  by  the  small 
square.  Thus  it  has  been  demonstrated  that  there  is  essentially 
perfect  agreement  between  the  analytically  predicted  eigenfrequencies 
and  those  determined  by  the  PETROS  3.5  code  provided  the  amplitude  of 
response  is  small. 

Turning  to  the  s = 2 family  of  vibratory  modes,  whose  lowest 
frequency  is  two  orders  of  magnitude  greater  than  the  fundamental 
flexural  frequency  in  the  unpressurized  case,  and  making  use  of 
Equations  (14)  and  (15)  it  may  be  deduced  that  this  is  an  extensional 
set  of  modes  which,  in  general,  also  entail  flexural  stresses.  The 
n = 0 mode  of  this  set  is  the  axisymmetric  breathing  mode  of  the 
cylinder.  The  s = 2 modes  exhibit  only  a slight  dependence  on  internal 
pressurization,  as  shown  in  Figure  17. 

The  importance  of  the  extensional  modes  in  computations  of 
accelerations  was  cited  previously  in  connection  with  Figures  7 and  8. 

To  better  define  the  contribution  of  these  modes  to  acceleration 
response  the  acceleration-time  data  for  6=0®  corresponding  to  the 
PETROS  3.5  run  represented  by  the  dashed  curve  in  Figure  5 was  analyzed 

by  use  of  the  Fast  Fourier  Trans form^^^^^  over  the  range 
80  ms  < t < 200  ms.  A portion  of  the  resulting  spectral  density 
function  is  plotted  in  Figure  18.  Since  the  range  selected  only 
slightly  exceeded  the  period  of  the  fundamental  inextensional  mode 
there  is  poor  resolution  of  the  amplitudes  of  the  first  few  s = 1 modes. 
Also  indicated  on  Figure  18  are  the  frequencies  of  the  rocket  motor 


W.  Cooley  and  J.  W.  Tukey,  "An  Algorithm  for  the  Machine 
Calculation  of  Complex  Fourier  Series",  Mathematical  Computations, 
vol.  19  (1965),  pp.  297-301. 

M.  Brenner,  "Three  FORTRAN  Programs  that  Perform  the  Cooley-Tukey 
Fourier  Transform",  Lincoln  Laboratory,  M.I.T.,  Technical  Note 
1967-2. 
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Figure  18,  Acceleration  Spectral  Density  at  0 = 0® 

calculated  by  use  of  Equation  (11)*.  These  results  appear  to  confirm 
that  the  (0,2)  and  (1,2)  extensional  modes  as  well  as  some  of  the 
inextensional  modes  are  definitely  present  in  the  acceleration  response 
of  Figure  7, 


Since  the  analytical  modeling  leading  to  this  equavion  did  not  provide 
for  any  circumferential  variation  of  density  an  "effective"  uniform 
density  was  employed  based  on  the  total  mass  of  motor  case  plus 
propellant.  This  approximation ^ which  "smears  out"  the  inhomogeneity 
associated  with  the  slots  in  the  propellant,  is  not  believed  to 
introduce  any  large  error  in  the  frequency  predictions  for  small 
values  of  n. 
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B.  Forced  Vibrations 


The  response  of  the  plane  strain  cylinder  to  the  prescribed  blast 
loading  has  been  analyzed  in  the  classical  manner,  representing  both 
the  applied  loading  and  the  radial  deflections  by  infinite  series  in 
the  free  vibration  eigenfunctions: 


n=0 

OD 

w(e.t)  cos  n0  (17) 

n=0 

The  coefficients  for  the  load  series  may  be  determined  from 

IT 

Po(t)  = 7 y*  Pg(0,t)  de  (18) 

0 

TT 

p (t)  = — f p (0,t)  cos  n0  d0  for  n > 1 (19) 

n TT  j e 

0 

When  the  specific  load  function  of  Equation  (1)  is  introduced  into 
Equations  (18)  and  (19)  and  the  integrals  are  evaluated  (numerically) 
the  values  of  Pj^(t)  shown  in  Figure  19  are  obtained.  The  intent  in 

presenting  this  Figure  is  not  to  present  quantitative  values  of  the 
coefficients  but  merely  to  show  that  their  variation  is  reasonable  and 
that  their  amplitudes  decrease  monotonically  after  full  envelopment  of 
the  cylinder  by  the  blast  wave,  indicating  a reasonable  rate  of 
convergence  for  the  series  of  Equation  (16) . 

As  a check  upon  the  adequacy  of  the  Fourier  representation  of  the 
applied  loading  given  by  Equation  (16),  plots  were  made  comparing  the 
initial  blast  loading  of  Equation  (1)  with  10-  and  20-  term  Fourier 
representations  at  various  times  during  and  after  the  diffraction 
phase  of  loading.  Three  of  these  plots  are  shown  in  Figure  20.  It 
may  be  seen  that  at  the  early  stages  of  envelopment  the  maximum  number 
of  terms  in  the  series  representation  is  desirable  whereas,  at  later 
times,  fewer  terms  would  be  adequate.  However,  for  results  to  be 
presented  subsequently,  the  full  20  - term  series  was  employed. 
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Figure  19.  Transient  Values  of  External  Load  Coefficients 

When  the  expressions  for  pe  and  w given  by  Equations  (16)  and  (17) 

are  substituted  into  Equation  (9)  to  effect  a separation  of  variables 

the  following  ordinary  differential  equations  for  the  coefficients 

w (t)  are  obtained: 
n 

“n  * - k)-  rj""  " ^]”n 

+ 3^k(l+q)n2|n‘*-  (2  - n^  + 1 - ^ | w^  = - aBk|p^  + 3(l+q)n2p^|  (20)  * 

n = 0,1, 2, 3, 


* ..... 

i4  dot  over  a variable  indicates  differentiation  wrth  respect  to  ttme. 
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Figure  20.  Fourier  Representations  of  Pressure  Distributions 
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psi 


These  equations  were  solved  in  the  conventional  manner  employing  the 
Laplace  transform  technique,  imposing  the  initial  conditions  of  zero 
displacement  and  velocity  at  t = 0 for  both  radial  and  tangential 
components  of  displacement.  Since  the  tangential  displacement  v was 
eliminated  in  the  derivation  of  Equation  (9)  (at  the  expense  of  raising 
the  order  of  the  differential  equation  for  w)  it  was  necessary  to  have 
recourse  to  Equations  (5)  and  (6)  to  establish  initial  values  of  w and 

w which  would  ensure  that  v(e,t)  = 0 and  4^(0, t)  = 0,  When  the  solutions 

n ot 

for  the  individual  are  inserted  in  Equation  (17)  the  total  radial 

response  component  is  obtained: 


t 

^ 5^  y “0,; 

- ’ 0 


sin  (Jjj  2“^ 


00  X 

ECOS  n0  /*  r*  'v 

~2 2 — / 

tvs'?  •' 


B(l+q)n2-u)2^^ 


Sin  w 

n,l 


sin  0)  -T  1 dx 
n,2 


g(l^q) 


pK.2 


cos  0 / Pj(t-T)  T dx 


The  last  term  in  this  equation  is  a non-vibratory  term  which  represents 
the  rigid  body  translation  of  the  rocket  motor  section  in  the  direction 
of  the  0=0  axis  produced  by  the  transverse  blast  loading.  It  may  be 
omitted  if  one  is  only  concerned  with  the  deformation  of  the  section 
relative  to  its  moving  mass  center.  The  circular  frequencies  appearing 
in  Equation  (21)  are  the  free  vibration  eigenfrequencies  defined  by 
Equation  (11).  If  the  tangential  displacements  v(0,t)  are  also  desired 
they  may  be  readily  obtained  by  substituting  the  result  for  w(0,t)  from 
Equation  (21)  into  Equation  (8)  and  performing  a single  integration 
with  respect  to  0. 
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A digital  computer  program  has  been  written  to  facilitate  the  use 
of  Equation  (21).  This  program  evaluates  the  convolution  integrals 
appearing  therein  by  numerical  integration  for  prescribed  pressure 
coefficients  p^^Ct),  such  as  those  displayed  in  Figure  19.  This 

formulation  is  quite  satisfactory  provided  w(e,t)  is  to  be  evaluated 
for  only  a few  values  of  time.  However,  if  w is  to  be  calculated 
for  numerous  values  of  t (e.g.,  to  obtain  data  for  a time-wise  plot  of 
w)  the  evaluation  of  the  convolution  integrals  over  the  range  0 < t < t 
becomes  costly  and  wasteful.  This  difficulty  can  be  circumvented  by  use 
of  the  following  computational  algorithm.  Equation  (21)  contains 
numerous  integrals  of  the  type 


/p  (t-T)  sin  0)  T dr  = / p (x)  sin  u (t-x)  dx 
^n^  ^ n,s  J *^n^  ^ n,s^ 


= sin  w t / p (x)  cos  to  x dx 
r\,s  J *^n'^  n,s 


- cos  to 


' t / p (x)  sin  to  X dx 
n.s  J ‘^n'-  n,s 


Let  the  time  interval  0 < x < t be  divided  into  increments  At  such  that 
mAt  = t and  introduce  the  notation 


/ 

(j-l)At 


(x)  cos  (0  X dx 
' ^ n,s 


J 

. f Pn'’’ 

n,s  J 


sin  to  X dx 
n,s 


IIIU  c 

= / Pn 


(X)  cos  to^^gX  dx  = 


s (t)  = f P (t)  sin  0)  T dx  = s"* 

J n n,s  n,s 

0 


Equation  (25)  may  be  written 


“n.s  ■ J 


(m-l)At 


Pn(T)  cos  * J “"n.s^ 


(m-l)At 


C™"^  + AC™ 
n,s  n,s 


= AC^  ^ + AC™  = y''  AC^ 

n,s  n,s  X V n,s 

j=l  j=l 

Thus  C (t)  can  be  evaluated  by  accumulating  the  sum  of  the  AC^  to 

° j,,S 

obtain  C™  and  then  only  having  to  calculate  AC™  (by  numerical 
"»s  n, s ^ 

integration  unless  the  interval  At  i^  taken  so  small  that  a cruder 

approximation  suffices)  for  the  current  At.  The  same  procedure  can  be 

applied  for  calculation  of  S .(t).  Then  Equation  (21)  can  be  rewritten 


The  computer  program  which  was  developed  for  evaluation  of  Equation  (21) 
has  been  modified  to  incorporate  this  algorithm. 


C.  Transient  Response  Correlations 

A comparison  of  deflections  at  6 = 0 obtained  with  the  PETROS  3.5 
code  and  with  Equation  (28)  of  the  analytical  solution  is  shown  in 
Figure  21  for  the  standard  blast  load  applied  to  the  bare  motor  case 
^ with  Pj^  = 6.89  MPa  (1000  psi) . The  blast  wave  arrives  at  the  pre- 

pressurized shell  at  t = 0,  driving  the  shell  at  6 = 0 inwards  and 
exciting  the  vibratory  response  which  follows.  The  agreement  between 
these  two  solutions  is  entirely  satisfactory  for  engineering  purposes 
in  this  case  of  small  amplitude  response  where  the  linear  elastic 
analysis  should  be  valid.  In  fact,  the  analytical  solution  should  give 
accurate  results  for  all  values  of  p.  to  the  right  of  the  "knee"  of  the 
curve  in  Figure  12.  Conversely,  one^would  not  expect  it  to  provide 
1 accurate  results  where  the  amplitude  of  response  is  large  and  this  is 

i borne  out  by  the  comparison  of  responses  for  p^  = 0 shown  in  Figure  22*. 


Figure  21.  Comparison  of  Predicted  Responses  for  Pressurized  Motor  Case 

For  this  comparison  (unlike  the  result  shorn,  in  Figure  5)  the  yield  stress 
of  the  motor  case  material  was  made  very  large  so  that  the  PETROS  2.5 
solution  remained  elastic,  thus  providing  a proper  basis  for  comparison 
with  the  series  solution. 


hi 
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Figure  22,  Comparison  of  Predicted  Responses  of 
Unpressurized  Motor  Case 

There  is  an  interesting  explanation  for  the  deviation  illustrated  in 
this  figure.  Note  that  whereas  the  usual  geometric  nonlinear  effect 
is  a "hardening"  phenomenon,  here  the  behavior  is  an  apparent  "softening" 
of  the  cylinder  when  the  amplitude  of  the  elastic  response  becomes 
large;  i.e.,  the  true  response  amplitude,  represented  by  the  PETROS  3.5 
solution,  has  a significantly  greater  value  than  that  of  the  linear 
series  analysis.  The  cause  for  this  behavior  can  be  related  to  the 
mean  component  of  the  external  pressure  shown  in  Figure  15.  Whereas 
the  PETROS  3.5  code  correctly  takes  account  of  the  transient  membrane 
stresses  induced  by  this  portion  of  the  external  load  the  analytical 
model  can  only  represent  the  effect  of  a static  internal  pressurization. 
To  show  that  this  is  the  correct  explanation  for  the  deviation  between 
the  two  curves  labeled  p^  = 0 in  Figure  22,  the  series  solution  was 

re-run  using  Pj  = ~ 0.115  MPa  (-  16.7  psi) . The  resulting  curve  (so 

labeled)  may  be  seen  to  be  in  much  better  agreement  with  the  PETROS  3.5 
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•r^ 


calculation*.  Of  course,  in  the  case  of  laige  internal  pressurization 
neglect  of  the  mean  external  pressure  is  of  no  consequence. 

The  effect  of  varying  the  blast  wave  intensity  on  deflection 
response  has  been  studied  for  the  bare  motor  case  with  no  internal 
pressurization.  The  results  for  the  maximum  deflection  (occurring  at 
0=0)  versus  the  peak  reflected  pressure**  are  presented  in  Figure  23. 
The  solid  curve  is  based  on  PETROS  3.5  calculations  in  which  the  motor 
case  material  was  required  to  remain  elastic.  The  analogous  response 
for  the  material  of  Figure  2 is  shown  by  the  short -dashed  curve,  the 
point  surrounded  by  a triangle  corresponding  to  the  "without  propellant" 
curve  of  Figure  5.  The  point  surrounded  by  a square  represents  the 
series  solution  for  p^  = 207  kPa  (30  psi)  and  p^  = 0.  It  may  be  seen 

that  the  solid  curve  is  tangent  to  the  straight  line***  connecting  this 
point  to  the  origin  for  small  p . Of  course,  as  noted  in  the 

max 

previous  paragraph,  the  series  solution  can  be  modified  to  give  a 
maximum  deflection  in  close  agreement  with  the  elastic  PETROS  3.5 
result  by  use  of  an  "effective"  p^  equal  to  the  negative  of  the  mean 

external  pressure.  However,  the  more  rigorous  PETROS  3.5  calculation 
is  to  be  preferred. 


VI,  CONCLUDING  REMARKS 

The  foregoing  analysis  of  solid  propellant  rocket  motor  response 
to  air  blast  loading  was  restricted  to  the  plane  strain  case  partly  to 
conform  to  present  design  practice  and  also  because  it  was  desired  to 
proceed  in  an  evolutionary  manner.  While  this  restriction  limits  the 


This  discussion  is  only  intended  to  identify  the  source  of  the 
dominant  nonlinear  effect.  For  satisfactory  quantitative  results 
the  PETROS  3.5  code  should  he  used. 

A A 

The  blast  wave  cannot  be  defined  by  a single  parameter.  The 
parameters  p , U,  ^ were  determined  for  a selected  p to  correspond 

to  the  blast  from  a one  kiloton  nuclear  device, 

AAA 

It  is  not  contended  that  the  series  solution  predicts  that  the 
maximum  deflection  is  a strictly  linear  function  of  p^  . However ^ 

^max 

since  the  peak  deflection  occurs  at  times  < < 2/X  and  X varies 
slowly  with  p the  response  predicted  by  the  series  solution 
^max 

should  be  nearly  linear  in  p^ 

^max 
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applicability  of  its  results  to  long  slender  motors  the  failure 
mechanisms  revealed  will  certainly  have  to  be  considered  when  more 
general,  three-dimensional  analyses  are  performed. 

The  finite  amplitude  elastoplastic  shell  response  code  PETROS  3.5, 
while  not  capable  of  fully  representing  the  response  of  a rocket  motor 
containing  propellant,  has  provided  useful  baseline  data  for  evaluating 
predictions  derived  with  other  models.  The  linearized  analytical  model 
provides  insight  into  the  meaning  of  the  results  of  the  nonlinear 
numerical  analysis  even  where  certain  assumptions  made  in  its  derivation 
are  violated;  e.g.,  the  large  amplitude  response  of  the  unpressurized 
motor  when  the  blast  loading  approaches  the  critical  buckling  load.  In 
addition  to  providing  an  independent  check  on  the  validity  of  the  PETROS 
3.5  solution,  the  modal  series  solution  is  preferred  for  long  term 
response  both  for  economy  and  because  it  does  not  entail  the  progressive 
discretization  error  of  the  time-marching  numerical  analysis.  Of  course, 
it  is  appreciated  that  neither  of  these  approaches  models  the  dissipative 
mechanisms  which  are  present  in  physical  rocket  motor  structures  but 
they  can  certainly  provide  an  upper  bound  to  the  response  parameters. 

The  response  of  a typical  rocket  motor  configuration  was  calculated 
for  the  limiting  situations  of  the  bare  motor  case  and  of  the  motor 
case  containing  the  complete  propellant  grain,  each  with  no  internal 
pressurization  and  with  the  pressurization  resulting  from  propellant 
combustion.  It  was  found  that  the  unpressurized  motors  had  a much 
larger  deformation  than  the  pressurized  motors.  Also,  for  the  motors 
containing  propellant  it  was  possible  to  obtain  an  estimate  of  the 
tensile  stress  to  which  the  propellant/case  bond  would  be  subjected. 

The  most  significant  result  is  the  quantification  of  the  greater 
vulnerability  of  rocket  motors  prior  to  their  ignition,  whether  on  the 
launcher  or  as  upper  stages  of  in-flight  missiles.  Consideration 
should  be  given  to  the  possibility  of  enhancing  the  survivability  of 
rocket  motors  when  introduced  to  a potential  nuclear  battlefield  area 
by  providing  a pre-pressurization  sufficient  to  avoid  the  large 
amplitude  response  to  air  blast  characteristic  of  unpressurized  motors. 

The  combination  of  the  PETROS  3.5  and  modal  series  methodologies 
provides  the  designer  with  a useful  approach  for  assessing  the  structural 
integrity  of  proposed  or  existing  rocket  motor  configurations  when 
exposed  to  a hostile  load  environment.  Several  extensions  of  the  work 
reported  herein  are  clearly  desirable.  The  finite  element  modeling  of 
the  viscoelastic  propellant  grain  presently  under  development  will 
provide  the  designer  with  required  information  concerning  critical 
stresses  both  within  the  grain  and  at  the  propellant-case  interface. 

A three-dimensional  analysis  model  is  also  needed,  especially  for  the 
shorter  rocket  motors.  For  an  adequate  vulnerability  assessment  of 
rocket  propulsion  systems  in  the  presence  of  nuclear  detonations  an 
analysis  including  the  nuclear  thermal  flux  deposition  properly  time- 
phased  with  the  air  blast  arrival  must  be  considered. 
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LIST  OF  SYMBOLS 
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Integral  defined  by  Equation  (25) 

Young's  modulus 

Circumferential  bending  moment/unit  length 

Outside  radius  of  motor  case 

Integral  defined  by  Equation  (26) 

Shock  front  velocity 
Mid-surface  radius  of  motor  case 
Sonic  velocity  at  ambient  conditions 

Natural  frequencies  of  cylinder 

Thickness  of  motor  case 
Summation  index 

(h/a)2/12 

t/At 

Circumferential  wave  number 
External  (blast)  pressure 

Internal  pressure 
Ambient  pressure 

Coefficient  of  n^h  term  of  series  for  p 

Reflected  peak  overpressure 

Side-on  overpressure  of  incident  shock 

(l-v2)p^a/(Eh) 

Mode  family  identifier 
Time 

Time  of  arrival  of  shock  front 
Circumferential  component  of  displacement 
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LIST  OF  SYMBOLS  (Continued) 


Radial  component  of  displacement 
Amplitude  of  n^^  term  in  series  for  w 

Integral  defined  by  Equation  (23) 

Integral  defined  by  Equation  (24) 

Time  increment 

12(l-v^)a‘* 

Eh  3 

Eh 

(l-v^)pha^ 

Circumferential  membrane  strain 

Circumferential  angular  coordinate 
Blast  wave  decay  coefficient 
Poisson's  ratio 
Mass  density 
Yield  stress 

Dummy  time  variable  in  convolution  integrals 
Period  of  the  (n,s)  mode 

Circular  frequency  of  (n,s)  mode 
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